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ABSTRACT 


Numerical integration and Monte Carlo techniques are 
used in the development of several models in order to 
determine the effect on probability of random detection 
of a merchant ship using speeds up to 90 knots by a 10 
knot submarine patrolling a back-and-forth barrier. A 
definite range law for detection is assumed. Individual 
encounter models are developed for ship tracks which cross 
the midpoint of the submarine patrol line at various angles. 
The models are extended to include the assumption of a 
normal-distribution of crossing points. Computer programs 
of the models, written in the FORTRAN IV language, are in- 
cluded. The results are applied in a numerical example. 

It is concluded that while increases in ship speeds do 
result in a substantial decrease in probability of detection 
by a submarine in the case of a single barrier transit, a 
speed advantage alone when applied to a typical transit of 
the North Atlantic will not appreciably decrease the overall 


detection probability. 
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I. INTRODUCTION 


A. BACKGROUND 
Pe Hace eee ol ome velopment 

New principles and concepts in ship hull design 
have made possible the development of ships capable of 
speeds heretofore considered impossible. In particular, 
the advent of surface-effect ship (air cushion and captured 
bubble), hydrofoil, and catamaran hull designs, coupled 
with advances in nuclear power, gas turbine, and water-jet 
propulsion systems, will have revolutionary effects on 
naval strategy and tactics, and world shipping patterns. 
That serious consideration to high-speed merchantmen is 
being made is evidenced by the fact that the U.S. Navy and 
the Maritime Administration share funding in the Joint 
Surface-Effect Ship Program which was established in 1966 
to determine the feasibility of 4,000 to 5,000-ton surface- 
effect ships capable of speeds in excess of 80 knots. Bell 
Aerosystems, Inc., and Aerojet-General Corp. have each 
received $1.5 million as a start on a 3-year program to 
build an experimental 100-ton surface-effect ship [1]. The 
fact that high-speed merchant ships which can enjoy large 
speed advantages over submarines are indeed in the offing 
indicates a need for reexamination of the abilities of 
future submarines to detect and kill these ships. fThis 
thesis makes a start by examining the effect of ship- 


Submarine speed ratios on detection probability 


is 


resulting from random encounters by submarines on fixed 
patrols. 
a. Literature Search 

As might be expected, the literature dealing 
with high-speed ships is relatively sparse; those papers 
that do make reference to ships of the future concern 
themselves with the hydrofoil small-ship types [2], or with 
the so-called super merchant ships which have such features 
as large (500,000-ton) displacements, and underwater 
tankers of speeds up to 40 knots [3]. Various Navy studies 
in the past which have attended to the question of ship 
speeds versus vulnerability have been limited in scope to 
surface ship speed ranges up to only 35 knots [4, 5]. One 


study in particular concluded that "...When ships are sailed 
independently... (and) when cost is considered there is only 
a Marginal gain in going from 20 to 24 Knots...." /1/ ee 
cost consideration at that time (1955) involved studies of 
conventional hulls. An earlier study [6] showed a decrease 
in percentage of independent merchant ships sunk from 10-12 
per cent when ship speeds were approximately equal to 
submerged submarine speeds (4-8 knots) to 1-2 per cent when 
ship speeds were double that of the U-boat submerged speeds. 
When hull-mounted sonars on merchant ships were considered 
[7], maximum optimal sustained speeds for independent 
vessels were in the range of 20 to 24 knots. 


In summary, neither the search of available 


studies [8], operations research papers, nor computer 
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searches for documentation [9] at the Defense Documentation 
Center, Washington, D.C., and at the U.S. Naval Postgraduate 
School, Monterey, California have revealed any literature 
regarding high-speed ships, and in particular the effect of 
merchant ship speed increases beyond 35 knots on vulner- 
ability to submarines. Indeed, all facets associated with 
these high-speed ships should be examined, e.g., noise 
radiation in water, the submarine fire control problem, 
effect of torpedoes, etcetera. 

2. Vulnerability of Convoys and Independent Ships 

As was clearly evidenced during World Wars I and II, 

convoys and independent merchants were quite susceptible to 
attack by submarines. Allied shipping losses from the 
German U-boat campaign during the first five years of World 
War II provide sufficient evidence of the extreme vulner- 
polity of merchant shipping [10, ll]. Considering che 
improvements in both submarine and submarine weapon system 
Capabilities since this era, the shipping losses previously 
experienced can be regarded as the minimum losses that may 
be expected ina future war [10]. Furthermore, it is a 
foregone conclusion that in the opening stages of a future 
general war escort numbers will be insufficient to provide 
adequate convoy protection for all merchant shipping [12, 
13]. It can be assumed, then, that the initial phases of 
a future war will see an aggressive enemy submarine force 
presented with many unescorted targets which can be attacked 


WLeh limeunity [10]. 
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To counter the possibility of sustaining heavy 
merchant losses at the outset of a future war, passive 
defense measures such as installation of sonar on merchant 
ships, routing to minimize losses, convoying, and increas- 
ing merchant ship speeds have been considered [{12]. This 
thesis considers the passive defense measure of increased 
merchant ship speed. 

cee bt Urem Coed enc mommoUloman ines 

The Type XXI submarine produced by the Germans in 
1944 influenced basic submarine design until the appearance 
of the nuclear-powered submarine in 1954 [14]. Since the 
introduction of the nuclear-powered submarine, great tech- 
nological strides have been realized in the areas of sub- 
marine hull construction, propulsion, sonar, torpedo 
weapons systems, and fire control systems. These advances 
have contributed significantly to the threat to merchant 
shipping posed by submarine forces of potential enemies. 
The modern nuclear-powered submarine can equal or exceed 
the speed of its surface-ship target or pursuer, and can 
remain submerged for days. Together with this endurance 
capability, the nuclear-powered submarine is able to main- 
tain a high sustained submerged speed in a relatively 


quiet environment. 4 .The limiting lines” of approachWee 
convoys and task forces steadily widen with increasing 
Submarine speed until there is no sector safe from attack 


bye awsOmeknet submarine. .@e" /2/ This quote was taken from 


a publication dated April 30, 1962; since this time a great 
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deal of research and development in the area of surface- 
effect ships has taken place. The Joint Surface-Effect 
Ship Program was initiated to consider large surface- 
effect ships capable of high speeds. It is conceivable 
that in the forseeable future surface ships with speeds 
far exceeding the capabilities of the nuclear-powered 
Submarine will take their place in the surface ship com- 
munity. Submarine speeds are not expected to exceed their 
present maximum for some time to come, thus permitting 
the high-speed merchant ship to enjoy a comfortable ad- 
vantage in speed. 
eee city er noe peleenies 

Heretofore convoys have been recommended as the 
best alternative for moving merchant vessels across the 
oceans when a submarine threat exists [15]. This is due 
primarily to the fact that for slow moving (less than 15 
knots) merchants, safety lies in numbers. The measure of 
effectiveness for analysis of merchant shipping has been 
percentage lost of the total number of ships in the 
convoy [16]. Studies completed during World War II indi- 
cated that approximately the same number of ships were lost, 
per attack, from large convoys as from small convoys. Hence 
the conclusion was to not only continue convoying, but to 
increase the size of the convoys as well. 

With the advent of high-speed surface-effect ships, 
the problem of independent sailing versus convoying should 


again be considered. In a study conducted by the 


yy 


Antisubmarine Operations Research Group Tenth Fleet in 
December 1943, the advantages and disadvantages that might 
result from independent sailings were discussed: "....Fast 
merchant vessels, with the exception of those carrying high 
priority cargo, have been able to continue sailing inde- 
pendently along the coast (Atlantic) much of the time." /3/ 
Here is a direct reference to a merchant ship speed advan- 
tage over the submarine. 

Other advantages of independent sailings are of 
two kinds: first, the independent merchant can operate more 
efficiently and hence deliver more cargo in a given time; 
second, the reduced requirements for escorts releases these 
surface, air, and submerged craft for offensive roles [ll]. 

5. Probabilities of Encounter 

The problem of random encounter between two units, 
sometimes referred to as the searcher and the target [l7, 
18, 19, 20], is examined in most of the classical operations 
research literature, Operations Evaluation Group reports 
and studies (which constitute a large part of the so-called 
"Classical literature"), and more recently by search theor- 
ists and game theorists. 

Probabilities of encounters, or detections resulting 
from applications of search plans, are discussed in much of 
the literature, ranging from Kimball and Morse [16] and 
Koopman [18], through the summary reports of the Operations 
Evaluation Group [14] and applications studies [2] to the 


modern search theorists, such as Pollock [19, 20]. 
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Data concerning encounters versus ship position 
displacement from reported submarine contacts [21] and 
frequency of encounters versus submarine position predic- 
tions [22] have been gathered and analyzed during the 
Second World War. Applications of game theory for the 
selection of convoy routes which would minimize ship- 
Submarine encounters were proposed in a study by Danskin 
of the Operations Evaluation Group in 1954 {15], which 
later appeared aS a paper in Operations Research [23]. The 
Random Encounter with no involved search effort has been 
studied by several authors; Koopman [{17, 18] considered the 
random encounter of an observer progressing at a constant 
velocity within a uniform distribution of targets of con- 
stant speed, as well as within a circular-normal distribu- 
tion of targets. The problem of random encounter between 
a surface force progressing at constant velocity with 
another unit travelling at a fixed speed, the position and 
heading of which are independently and normally distributed, 
has been examined by Dobbie in 1945 [24]. 

Literature concerning the effect of speed on the 
random encounter is limited almost exclusively to government 
studies and reports. Kittel, in 1944, analyzed data con- 
cerning sinkings and sightings by submarines to give an 
early indication of the effect of speed on the safety of 
ships [6]. The results were discussed by Sternhell and 
Thorndike in 1946 [25], and by Winston in 1955 [4]. These 


results, in turn, were summarized by the Operations 


ig 


Evaluation Group in 1957 [12]. More recently, the general 
problem of the effect of speed on the vulnerability of 
independent merchant vessels was examined by Neufer in 
1961 [4]. Indeed, Neufer appears to have been the first 
to attempt to determine the explicit survival probability 
of a merchant vessel in transit from one port to another as 
av runctven of 1ts speed. 

Einaldivjuethnemprebabinity oihvencounterubyanunikteon 
a barrier patrol has been examined by at least two authors. 
Koopman, in his Search and Screening, examined the problem 
using the underlying assumption of a uniform distribution 
of points of intersection of the target's track with the 
back-and-forth barrier patrol line. /4/ In particular, 
Koopman showed that, with a definite range law, the 
probability of detection by the back-and-forth barrier 


patrol is given by 


a in; alld 
Pe 1) - ve L/L), 


ae ome 1 see 
| 
a foreore 1/ (27 eo 
where 


A = ratio of the barrier patrol line length to 
sweep width 
6 = ratio of target speed to barrier unit patrol 


speed 
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Further, Koopman showed that, for all values of A, the 
back-and-forth barrier patrol is preferred to the crossover 
STElgure 61igitejaey pe Dadkrol for all valves of ¢ 2.17% /5/ 
Van Train [26] used a war gaming method in a limited way to 
determine the probability of detection of a transitor 
intercepting a multi-unit, single or multi-line fixed 
barrier, in which each submarine was constrained 

1) to remain in the center of an assigned zone, and 

2) to conduct intercepts using headings parallel to the 
barrier line. The submarine's course, and position on a 
"starting line" were assumed to have been normally- 


Gilcteibuted. 


hee LAE PROBLEM 
The basic encounter problem is that of determining the 
probability of detection of a surface ship by a submarine 
on an assigned back-and-forth barrier patrol, such as that 
discussed by Koopman [18]. The following assumptions were 
made for the models: 
1. Assumptions for the Basic Random Encounter Model 
a. The submarine patrols a geographically-fixed 
back-and-forth barrier centered on a known surface ship 
Eeack Or COnVvoy route. 
b. The submarine patrols at a constant speed, v. 
c. A definite range law for detection is used. 
d. The submarine patrol length, L, is four times 


the detection radius. 


Zu 


e. The submarine position at any time, T, is from 
a Undsormedis tribution. 

fe The point Of Vimtercscection of the Shipeeseaack 
and the patrol line (or extension of the patrol line) is 
from a normal distribution, with standard deviation of 1/6th 


of the length of the submarine patrol line. 


Caw METHODS OF SOLUTION 

Two methods of solution to the models were devised; the 
first was a numerical integration technigue, i.e., simula- 
tion by systematic sampling. The submarine's initial 
position was successively incremented by 500 yards (0.25 
nautical mile) and the minimum lateral range (CPA) between 
the submarine and the surface vessel was computed by 
numerical vector analysis. The baSic model considered 960 
initial submarine positions for each speed ratio considered; 
twelve ship-submarine speed ratios, ranging from 1 to 9 were 
examined. The basic models also examined nine angles of 
intersection of tracks, ranging from 90 down to 30 degrees. 
The extended model considered ten specific points of inter- 
section and eight speed ratios. 

The second method may be classed as a Monte Carlo method, 
at least within Brown's definition: "...any procedure which 
involves the use of statistical sampling techniques to 
approximate the solution of a ... physical problem...." /6/ 
The Monte Carlo method utilized virtually the same computer 
program subroutines as the numerical integration approach, 


but selected the submarine's initial positions from a 


Zz 


uniform pseudorandom number (URN) generator [27], and 
selected track intersection points using the Box-Mueller 
method of normal random number generation [28]. For this 
technique, the basic model took 1,000 samples for each of 
nine speed ratios per crossing angle, and for 17 crossing 
angles. The extended model took 10,000 samples for each 
of nine speed ratios. The total number of vector 
("maneuvering-board") solutions, then, exceeded 500,000; 
the usefulness of the computer in this analysis is obvious. 
1. The Random Encounter Model 

Consider a submarine patrolling a barrier of length 
L miles at a constant speed, v. Assume that the position 
of the submarine on its patrol line at any time, T, iS a 


uniformly distributed random variable, Y Without loss of 


5° 
generality, then, the problem can be considered to begin at 
time t when a surface vessel on a course a, projected to 
intersect the Submarine patrol line at its center, is a 
distance x / COs a from the point of intersection. Define 
the distance a to be equal to the range of the submarine 
detection (assuming a definite range law), R, which in the 


initial case was considered to have been 30 miles. Refer 


to Figure 1 and the definitions summarized below: 


> E in ; 

u=u cos @1+u Sin a Jj = ship velocity; 

> o Stee 

v=etv gj = submarine velocities; 

> > > : : 

W, = vV,- u = relative velocity of the 

1 1 ; 

submarine when it 1s on the 
upleg; 


o3 


Submarine 
Patrol 
Line 


(Xx, Vo 


p’"p 





Figure 1. Vector Diagram of the 
Fncounter (Relative Space) 
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W. = Va- U = relative velocity of the 
Submarine when it is on the 
downleg; 


| 
x 
i 


L - Y/Y = time remaining for the 
submarine to be on its 
initial leg before course 
reversal; 


Gi: 
Il 


(Y = y yy, = time remaining for the sub- 
P marine to be on its initial 
downleg before course 
reversal; 


T = L/v = length of time required for 
the submarine to travel 
entire length of its patrol 
line; 


= x i+ YyJ = initial position vector of 
P the submarine from the ship; 


ON 
II 
Mo 
ae 
= 


fe = position of the submarine 
when tt makeserceseecourse 
reversal. 


Zee Oren, cOsleCeormime Prep aor sey, Obl Encounter. 
(1) Select the initial value of ty SUC that a Y 5° 


(2) Specify the initial direction of submarine velocity; 


Ceol lt V. 


(3) Compute the values of t, ti, Ri, and R; set w= Wy - 
(4) Compute the dot product Ry °W = Z: 
(A eo TE Fa onOe CPA = R.; Gow “step i5) > 


et 
(b)> LE 2 < ©, compute the det product Rotw 


Il 
N 


1) re cee ee — aR 
2) If z' > 0, the CPA lies between =F and Roe and 


must be determined; 
a) set a = 1/2; 
b) compute y = (aw) /|w|; 


—-> —> —> 
c) compute the dot product (Ry + y)-w= 2"; 
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1s) tt ezeeeer add 1/2 to a, andmgortce 

step 2)b) above; 
: 7g > > 

22) Zee Ehen ClA = Ry + Vee, EG 
step (5). 
(53) (a) ats | CPA| < R, a detection occurred; go to step 
(1) if the number of samples taken is greater than N, 
otherwise go to step (6). 


(2) Lt |CPA| < R, there was no detection on this leg: 


1 fo = es . 
yi Woe Wie set Ww = Wo; 
> ee eee mtg: 

) 1f we Wo, set W = wii 


3) go to step (4). 

(6) P(Encounter) = # detections/ N. 

The general flow diagrams describing the logic for the 
models are contained in Appendix B. Note that there are 
two models, namely, the Basic Model, in which the ship track 
crosses the midpoint of the submarine patrol line at various 
angles, and the Extended Model, in which the crossing angle 
is fixed at 90 degrees, but the crossing position is varied 
along the track (and its extension). Appendix C contains 
the actual computer instructions used to generate the data 


analyzed herein. 


Deo OBSERIPTION CEs tbe COMPUT Rao ys... EM 

The computer system used in the development of this 
thesis was the high-performance IBM System/360 Model 67, 
which is located at Ingersoll Hall, U.S. Naval Postgraduate 


school, Monterey, California. This sophisticated data 
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processing system employs a FORTRAN G compiler, along with 
the standard software and hardware associated with the 360 
SYStem. 
Programming was accomplished using the standard FORTRAN 
IV language as proposed by the American Standards Associa- 
eHon “S.4.3 EPORTRAN Working Group [29] . 
Approximate execution times for the four programs in- 
cluded in this thesis were as follows: 
Basic Model: 45 minutes; 
Extended Model: 30 minutes; 
Random Sampling, Midpoint: 110 minutes; 
Random Sampling, Normally-distributed crossing: 


46 minutes. 
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ll. THE RANDOMSENGEGUNTER -— ANALYSiee 
BY NUMERICAL INTEGRATION - 
i RANDOM ENCOUNTER MODEL (FIXED CROSSING POINT) 
Wa DESCrI DELON Oise e Ode! 
The algorithm described in section 1.B. was used, 


except that the values of Y and the initial direction of 


ae 
the submarine were determined using numerical integration 
techniques as follows: 

The submarine's initial position for the first 
sample was established at the bottom of the patrol line, 


with the initial direction "up." The relative motion prob- 
lem was solved numerically to determine if a detection 
occurred for this configuration. For each subsequent 
sample, the submarine's initial position was incremented 
by 500 yards (174th of a nautical mile), until each 500 
yard increment in the initial direction had been sampled, 
up to the initial position lying on the patrol line extrem- 
ity. The submarine direction was then reversed and the 
procedure continued as before except that the succeeding 
initial positions resulted from 500-yard increments in the 
new direction. The number of samples, then, was 8 x (length 
of the submarine patrol line in nautical miles). The basic 
parameters of the model from which most of the data in 
Appendices C and D were obtained were 

Submarine patrol line length: 120 miles 

Submarine patrol speed: 10 knots 


definite range law detection range: 30 miles 
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These parameters were used in both the basic and extended 
models for the simulation by numerical integration, and for 
the Monte Carlo simulations. 
2. Results 
The data resulting from the simulation by numerical 
integration are contained in Appendix A; these data have 


been plotted in Figure 2 and Figure 3. Figure 2a shows how 


en abet! ae Lt Mel Fin OL ETOP pate — 





the detection probability varied with the ship-submarine 
speed ratio for a ship crossing the patrol center at various 
angles. The points at which detection probability was 
reduced by 90% and 95% have been connected by the two curves 
labelled "90%" and "95%," respectively.” These curves are 
shown in Figure 2b. Figure 3 is an expansion in scale of 
Figure 2b, showing greater detail. Figure 4 displays a plot 
of the Probability of Detection vs. Track ar ery Angle for 
a midpoint crossing. The behavior of P(D) as a function of 
a is of little practical interest for the case where a 
submarine commander has freedom to select a patrol line 


orientation with respect to a convoy route; a patrol 


“The results plotted in Figure 2 seems to indicate that 
a barrier unit's best tactic would be to make the patrol line 
cross the ship's track at a shallow angle. This would 
indeed be the case if the ship's track were known with cer- 
tainty, aS was assumed in the basic model discussed here. 
The more realistic assumption of uncertainty in the know- 
ledge of the ship's track makes the choice of a shallow 
angle undesirable from the viewpoint of the barrier unit. 
The design of an "optimal" barrier patrol line, the crossing 
angle of which depends on, say, the probability distribu- 
tion of a convoy route displacement, appears to be an 
interesting problem, but is not considered herein. 
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Submarine Speed: 10 kts. 
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a = Crossing angle in degrees 
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Figure 4. Probability of Random Encounter vs. 


Track Crossing Angle (Midpoint Crossing) 
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perpendicular to the expected surface ship track would 
normally be selected. Figure 4 could be useful if the 
submarine were restricted for some reason to orient its 
patrol line at an angle other than 90 degrees to the convoy 
route. The effect of relaxing the assumption that the 
surface ship track crosses the patrol line at the midpoint 
of the latter is examined in the next section. The assump- 
tion of a submarine patrol line length of 120 miles (four 
times the detection range) agrees with Neufer's "total 


attack width" which is derived in his paper. /7/ 


B. EXTENSION OF THE MODEL TO VARIABLE CROSSING POINTS 
1. Description of the Model 


The basic model contains the expression 


= x tana - L/2 
Yn D / 


to insure that the ship's track will bisect the submarine's 
patrol line. It did not seem unreasonable to assume that a 
submarine normally attempts to center its patrol on the 
expected convoy route; that is, the patrol is centered on 
an expected value. A normal distribution of tracks, each 
perpendicular to the submarine patrol line, was assumed for 
the extension of the model, so that Y is a normally- 
distributed random variable. In terms of the coordinate 


system used herein, 


E(Y,,) = -L/2 
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The logic used in the computer program is illustrated in 
Appendix B. The results of the computer calculations are 
plotted in Figure 5. 
2 eeenak Sis or Results 

The flattening of the curves for the speed ratios 
op = 1 and p = 9 suggested the possibility of two simple 
limiting shapes, namely the shapes resulting when p = 0 
and p = » As p> 0, 1.e., when the ship's speed u is such 
that u<<v, it can be seen that the probability of detection 
will approach unity. More precisely, for a given detection 
radius, submarine speed, and submarine patrol line length, 
there exists some value of ship speed such that the sub- 
marine will always have an opportunity to close within 


detection range. Let 


L = length of the submarine patrol line 

Vy = y-coordinate of the bottom limit of the sub- 
marine patrol line 

u = ship speed 

V = submarine speed 

0 = ship-submarine speed ratio 

ie = detection radius (definite range law) 

P(D) = probability of detection 


Then, for all u, such that 


22 / sly vee; 
ee, for all p*, such that 


e/a; 


5 
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Figure 5. Probability of Detection vs. Crossing 
Position (Perpendicular Tracks) 
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lim P(D) = 1 


let < L t+ 
i pea s i 


For the case of p > ©, i.e., Submarine speed approaching 
zero, the probability of detection is simply the probability 
that the submarine's position (which has been assumed to be 
uniformly-distributed along the patrol line) falls within 


distance r of the ship's track. Refer to Figure 6 below. 


ioe ate ee at 






: P(D) = (x + Ay)/L ~Y 
: —- P(D) = (x - Ay)/L yo 


P(D) = 2r/L - = 


me aE Ge NR a Re ci te —_ = —_ 


er ee tal a ET SE ee te 
<< on “-—- - we + mem eee : ‘oh me >| 


Figure 6. Probability of Detection vs. Crossing 
Position for Zero Submarine Speed 


It can be seen from Figure 6 that if the ship's track 
crosses the patrol line at a distance from the patrol line 
end-points greater than or equal to r, then the detection 


probability is simply 


Daas 


oy. 


If the ship's track crosses the patrol line at a distance 
less than r, say Ay, from a patrol line end-point, then the 


detection probability is given by 


(aera) 7 L 


Finally, if the ship's track does not cross the patrol line, 
but passes within detection range r, say Ay, of one of its 


end-points, the detection probability is given by 


ieee 1G 


Using the convention that the coordinate system origin lies 
on the ship's position at the instant that the ship becomes 


"detectable," the above argument may be formulated as 





follows: 
faeps, Hae ek + r, -r) 
| xv - min(/y.|,|/L + y } san [(y,) (y, a) oy 
) L 
P (D) = | (i 
ypeti-r xr] [-1 - r, -L + rj} 
| 7 ypetir,o)U Cor ge I ral 
The limiting case result, as presented herein, appears at 
first glance to disagree with the limiting case result ob- 
tained by Koopman; Koopman's result for a zero-speed 
barrier unit is P(D) = 0.333... /8/, and the result from 
the geometric integration for pp = ™ is 0.488. This latter 


value is a consequence of the assumption of normality; if 


Koopman's underlying assumption, namely a uniform 


38 


distribution of crossing points, were to be applied to this 
model, the limiting probability of detection would be 
(area under curve for eP = ~)x(density for uniform distri- 


bMeion) or, 


CORE a8) (0.5) (60) _ 60 _ 
180 + 180 = ron OF sere 0 16 
which is identically Koopman's result. The model expressed 


by equations (1) is, therefore, consistent with Koopman. 
A simple FORTRAN computer program was written (Appendix F) 
to compute the detection probability using equations (1); 


the results are plotted in Figure 7. Figure 7 also contains 





a plot of the computed results for comparison with the 
limiting cases. The agreement with equations (1) is evident. 
The behavior of Probability of Detection vs. Crossing 
Position in the region of ship-submarine speed ratios 
greater than three is unexpected (one would expect the 
curves to remain uninodal). At values of OP such that 0>4, 
the probability of detection resulting from a patrol, the 
center of which was offset from the expected track crossing 
by about 1/6th of the length of the patrol line, resulted in 
a higher detection probability than when the patrol was 
centered on the track crossing point. There data indicate 
that a submarine may effectively reduce the speed ratio by 


properly positioning itself. 
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COMPARISON OF PEAKS ON PROBABILITY CURVES 
FOR SHIP-SUBMARINE SPEED RATIOS OF 5 AND 9 


SPEED P(D) FOR CENTERED P(D) GOR (OnE -CENTERED 
RATIO PATROL LINE PATROL LINE (20 MILES) 
5) ~5104 -5604 
9 so0 22 Pee 


Table I, extracted from the data, reveals that the prob- 
ability of detection of a ship with a nine-to-one speed 
advantage over the patrolling submarine by a submarine with 
a 20-mile offset patrol line is greater than the prob- 
ability of detection of a ship with a five-to-one speed 
advantage by the same submarine with its patrol line 
centered on the crossing track. The table shows that with 
a 20-mile offset, the submarine would enjoy a higher 
probability of detection than it would for the case of a 
target with p = 5 and no offset. Figure 8 employs a con- 
tracted horizontal scale which serves to emphasize the nodal 


behavior discussed above. The conditional probability, 
P (Detection |Y_, = y), 


was determined by numerical analysis, and is plotted in 
Figure 7 for six values of p. The track crossing point was 


assumed to be normally-distributed, that is, N(0,400). Now 
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Figure 8. Probability of Detection vs. Track Crossing 
Position (Perpendicular Tracks) 
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P(Detection) = P(R= rr), 


where re[0,~), and R is the random variable representing 
the distance between the submarine and the surface ship. In 
the problem at hand, r = 30 miles, and since a definite 


range law was assumed, a detection occurs when R = 30, hence 
P (Detection) = P(D) 
1s a mass function. From Baye's formula, 


ERiy Cie f(r] ¥,, Saat LUO (2 
Pp p 


Integrating over the range of y, 


co oc 


f ta,y (fevidy = fo f(r|Y¥, = y)fy (yay . (3) 
P P 

The right side of equation (3) is the expectation of the 

conditional density f(r|Y_ = y), and the left side is the 


p 
Marginal density, in this case the mass density, of r; hence; 


© 


P(D) = pp(r) = J f(r a = y)f, (y)dy . (4) 
em OO p 


Since the conditional density in equation (4) was derived 
numerically, and no explicit mathematical expression was 
determined, the integral indicated was performed by geo- 
metric integration technigues. Appendix E shows how the 
area under the curve 


Z= f(r] ¥,, = ee 
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was obtained. Equation (4) may be rewritten as 


CO L/2+r 
f £(r|¥ =y)fy Cy dyea2 £(x|¥ = y)f, (y)dy, (5) 
—o p O p 
since 
f(r[¥, = y) = 0, Y Yel= 0/2 =e 2), 


and the integrand is symmetric. 

Figure 9 is a plot of the function to be integrated; 
BeomclMes Of 9 wSed exe 0, lls, 32, 3, °5, 9, and ©. Dine 
results of the geometrical integration are listed below in 


Hale TT, and are plotted in Figure 10: 


TABLE II 


GEOMETRIC INTEGRATION RESULTS 


SHIP-SUBMARINE 


SPEED RATIO Je) {UB} 
0 noe 
I .9405 
he 5 79) 
2 oc LG 
3 Sol oul 
5 25233 
9 POS 
& ~4885 
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Cursory examination of Figure 10 reveals that for speed 
ratios less than four, the assumption that a fixed crossing 
point will yield more pessimistic results in terms of 
higher probabilities of detection than those arising from 
the assumption of normally-distributed cross points is not 
valid. The length of the submarine patrol line was 

varied from 120 miles down to 30 miles in order that the 
effect of changing the ratio of patrol line length to detec- 
Elon range might be examined. Hagure 11 shows the results 
for a patrol line length of 90 miles (three times the 
detection range). Note that the binodal behavior of the 
probability curves which was present when a patrol line 
length of 120 miles was used is absent in this case, except 
for the situation where o0 = 9, where the binodal behavior 
is incipient. The nodal phenomenon observed indicates a 


possible area of interest for additional research. 
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Figure ll. Probability of Detection for Various 
Speed Ratios vs. Track Crossing Position 
(Perpendicular Tracks) for L=90 miles 
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LIL Se THE RANDOM - ENCOUNTER = SEMULATION 
BY MONTE CARLO METHODS - 
A. MIDPOINT CROSSING MODEL 

The encounter model developed by deterministic methods 
was used to build the models and computer programs for the 
stochastic simulation. All assumptions pertaining to the 
deterministic model were included in the Monte Carlo ap- 
proach. 

The submarine's initial position on his patrol line was 
determined for each sample taken by selecting from a uniform 
(0,1) pseudorandom number generator. This procedure agreed 
with the assumption that the submarine was equally likely to 
be anywhere on his patrol line at any time T. 

The submarine's initial up or down heading on its patrol 
line was determined by selecting a pseudorandom number from 
a uniform (0,1) variate generator. If this number was less 
than 0.5 the submarine was assigned an initial up heading; 
if the number was greater than 0.5 the submarine was 
assigned an initial down heading. If the random variate was 
equal to 0.5 another uniform random variable was selected 
and the above procedure was repeated. 

The midpoint crossing of the patrol line by the merchant 
was arranged by positioning the bottom of the patrol line in 
such a manner as to insure that this phenomenon occurred. 
Section (11) gives the analytical explanation dealing with 


the development of the deterministic approach. 
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The results obtained using the Monte Carlo technique are 
presented in Appendix A. Graphs depicting the results can 
be found in Figure 12. It should be noted that the results 
obtained utilizing the stochastic simulation agree with 


those obtained from the deterministic approach. 


Bee TENSION TO NORMALE pts. teUrTED CROSSING POINTS 

The model and computer program for normally-distributed 
crossing points of the convoy route about the center of the 
submarine patrol line is an extension of the midpoint 
crossing model. 

The submarine's initial position and heading were 
determined as for the midpoint crossing model. The point 
of intersection of the submarine patrol line and the convoy 
route was determined by generating a normal N(0,400) 
pseudorandom variate using the Box-Mueller Method [28]. The 
mean of the normal variate corresponded to the center of the 
patrol line. The problem was then generated in a manner 
Similar to that of the midpoint crossing model. 

A relatively small sample size of 10 normal variates, 
one for each 1,000 uniformly-distributed initial starting 
positions for the submarine, was generated due to the time 
limitations on the computer system. This sample size took 
approximately 46 minutes to execute. The computer program 
which appears in Appendix D has been written using a larger 
sample size in an effort to aid in any future duplications 


@irthne Simulation. 


50 


(enbtuyoe3 OTrAeD sqUOW) STbuy aynoy FO AszZowerzeg yRIM 
oTZey peeds qns-dtys *sA ARTTIqeqord uoT_Zoe3eq “ZT eAanhbtYg 


A/n=d oOT}ZeY pseds 
6 8 L 9 G 7 € G 


o Gk =0 


oS P= 


o 0G=0 


0 GG=0 


oV9=0 





ATT tTqeqorg uot}o.9390q 


51 


The results obtained using this simulation technique 
are tabulated in Appendix A. A plot of the results appears 
in Figure 13 and Figure 14. The results and analysis of 
this situation are similar to those obtained from the method 
of analysis by numerical integration in the extension to 


variable crossing points as is evident from Figure 14. 
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IV. AN EXAMPLE OF A NORTH ATLANTIC CROSSING 


A. ASSUMPTIONS 

It was assumed that the major western North Atlantic 
seaports are New York, Norfolk, Windward Passage, and 
Puerto Rico (Panama Canal traffic); the major eastern 
North Atlantic seaports were assumed to be the Faero Islands 
(representing the straits between Iceland and Scotland), 
the English Channel, and Gibraltar. This assumption results 
in twelve possible Atlantic-crossing convoy routes. These 
possible convoy routes were usSed as a basis to determine 
the following assumptions: 

1) allocation of enemy submarines 

2) lengths of submarine barrier patrol lines 

3) angle of orientation of the submarine patrol lines 

to the convoy tracks. 
A consequence resulting from the assumption of specific 
convoy routes is a deviation from the assumptions made by 
Koopman [18], Neufer [5], Dobbie [24], and others which was 
that of a uniform distribution of submarines throughout an 
area. Figure 15 is a sketch which illustrates the rationale 
OG Bie meetin ois made herein. 

Given that the only possible convoy routes are as shown 
in the er aure: an enemy might be expected to allocate his 
anti-shipping submarine patrols as follows: 

1. One submarine to patrol each port (since ASW 


activity would be expected to be heaviest near these ports, 
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greater numbers of patrols at these locations was not 
assumed). 

2. A submarine patrol at each expected juncture of 
convoy routes (this effectively places a barrier on each 
of the convoy routes which merge at the junction). 

3. No segment of convoy track greater than 500 miles 
to be unopposed by a barrier. 

The preceding assumptions require that an enemy have 
the capability to maintain 70 submarines on station in 
the North Atlantic. Note that this distribution results 
in 105 barriers against the possible convoy routes assumed, 
even though only 70 submarines are employed. 

Another assumption made was that the submarine barriers 
would be oriented according to the following rules: 

1. Patrolling of straits to be oriented so as to 
minimize the patrol line length. 

2. Patrolling of seaport exits (entrances) to be per- 
pendicular to, and centered on, a bisector of the widest 
angle formed by two departure courses. 

3. Patrolling of route junctions to be along the 
bisector of the largest angle formed by the intersection 
of the two convoy routes, and to be centered on the 
expected junction. 

4, All other patrols to be centered on, and oriented 
at right angles to the expected route. 

Additionally, all submarine patrolling speeds were assumed 


to be 10 knots, and all patrol lines were assumed to be 


Dil 


120 miles in length (unless shortened by geographical 


limitations). 


Bey APPLICATION. OF THESMODELS 


1. Orientation of Barriers on the Norfolk-English 


Channel Route 

The route between Norfolk and the English Channel 
is approximately 3,500 miles long, and intersects the routes 
from New York to Gibraltar, Windward Passage to the Faero 
Islands, and Puerto Rico to the Faero Islands. The rules 
assumed required that the barrier patrols on this convoy 


route be given the following parameters: 


TABLE “diet 


PARAMETERS OF THE BARRIERS FOR THE EXAMPLE 


BARRIER 

+ BARRIER LOCATION ag L 
Norfolk 5 120 miles 
N.Y¥.-Gibraltar junction 12 120 

3,4 Along the route 0 120 

5 W.P.-Faero junction 10 120 

6 P.R.-Faero junction Hey 120 

i English Channel 
(Land's End-Ushant Is.) 0 90 


2. Solution to the Example Problem 


In general, 


P (Detection) = 1 - P(No Detection), 
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and in this case, 
P(No Detection) = P(No Detection by Barrier #1) 

xP (No Detection by Barrier #2) 

KX fore 

xP (No Detection by Barrier n) , 
assuming stochastic independence for simplicity. In this 
example, n= 7. Since the probability of detection is a 
function of the angle of track a, the patrol line length L, 
the detection range r, and the ship-submarine speed ratio 
0, the probability of detection for barrier i may be 


represented as 


P. (D) = P(D;a;,L;,,r;,0;) ’ 


and 


P(D) = 


5 [1 = P(D;a,,L,,r;,9;)] 


a 


lan 


According to the assumptions made, 


r, = Fee FYZ SOL 
= eo Boor ae See 
and, 
7 
P(Detection) = [I [l - P(D;a,,L,,r,p)} (6) 
i=l 


The probabilities of detection for each of the barriers 


listed in Table 3 are tabulated in Table IV. 
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Using the results from Table IV in equation (6), 
i= (| o-62 i ee Dl po," 1 20,0 = 0,5 = one 
we (Dest — 90rd Ommse— 0g 


The results are tabulated in Table 5. 
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Table V indicates, ceterus paribus, that speed advantage 
alone for the surface ship will not lessen the probability 
of detection significantly by patrolling submarines for a 


typical North Atlantic crossing. 
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V. SUMMARY AND CONCLUSIONS 


A. MONTE CARLO VERSUS NUMERICAL INTEGRATION METHODS 

Figure 14 shows that close agreement exists between the 
Monte Carlo and numerical integration techniques; the numer- 
ical integration method is more economical in terms of 


computer usage time. 


B. THE INDIVIDUAL RANDOM ENCOUNTER 
1. Midpoint Crossing at Various Angles 
Detection probability decreased from unity at P<l 
for all angles of intersection between the ship's track and 
the submarine barrier patrol line, to 0.5 at 9p > 9 (for 
intersection at right angles). Ninety-five per cent of the 
decrease occurred at speed ratios of eight and below for 
all angles, and for speed ratios of 4.5 and below for inter- 
sections at right angles. Stated another way, for assumed 
Submarine speed of 10 knots, and using detection probability 
as the sole criterion for selecting ship speed capabilities, 
merchant ship speeds greater than 40 knots do not seem 
warranted, 
2. Normally-Distributed Crossing Points 
a. The assumption of normally-distributed crossing 
points resulted in lower detection probabilities for ship- 
Submarine speed ratios less than four, than did the assump- 
tion of a fixed crossing point at the center of the patrol 
line. For speed ratios greater than four, both assumptions 


yielded essentially identical results. 
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b. If ship-submarine speed ratios equal to or 
greater than four are to be considered, then the binodal 
behavior of the probabilities of detection, illustrated in 
Figure 8, should be examined further. It appears that 
this behavior may also be a function of the ratio of sub- 


marine barrier patrol length to assumed detection range. 


C. NORTH ATLANTIC TRANSIT EXAMPLE 

When six or more submarine barriers are crossed by the 
merchant ship's track, the reduction in P(Detection) 
resulting from increasing the merchant ship's speed is 
negligible. It is evident that, from the viewpoint of 
decreasing random detection probability, a large speed 
advantage is not the answer. If a large speed advantage 
by a merchant ship over a patrolling submarine increases 
the survivability of that merchant ship, then the increase 
in survivability is due, not to decreased detections, but 
to other factors, e.g., reduced submarine weapon system 


effectiveness or reduced exposure time. 
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MONTE CARLO TECHNIQUE 
PROBABILITY OF RANDOM ENCOUNTER 
MID-POINT CROSSING 


DEFINITE RANGE LAW 


DETECTION RANGE: 30.0 MILES 


LENGTH OF SUB PATROL LINE: 120.0 MILES 
X=DISTANCE TO PATROL LINE: 30.0 MILES 


SUBMARINE SPEED: 10.0 KNOTS 


NO. OF SAMPLES TAKEN FOR EACH SCENARIO: 1000 

BVG. NO. SHIP SPEED RELATIVE SPEED ALFA PROBE 
LEGS SPEED RATIO UPLEG DOWNLEG DETECT He 
1.41 Oe ie ea. 1 ao 0.0 1.0000 
1.34 200 20 22.4 22.4 OF0 0.8200 
eS 3050 0 5.6 0 CEO 0.6270 
ile 40.0 4.0 Aas 2 41.2 0.0 0.5470 
reas DOR 550 a. 0 20 0.0 0.5 Zag 
eel 7 60.0 6.0 60.8 60.8 0.0 0.5060 
eel 70. 0 Uo 707 Obed 00 0.4810 
ee 80.0 8.0 80.6 80.6 0.0 0.4810 
ea 3 90.0 o.0 90.6 90.6 0.0 0.5070 
ena 2 yo) 1. eo ese 14.7 020873 1.0000 
er 2 207.9 2.0 21.6 23750 020873 0. 7650 
a2 5 30.0 3.0 30.8 32.4 OF08 73 0.6420 
9 40.0 4.0 40.4 42.1 0.0873 0.5460 
mee 30-10 2.0 a0 1 Sao 0.0873 0.5220 
Lo 60.0 6.0 60.0 61.7 OR0873 0.5020 
eS 70.0 7.0 Soo eG 0.0873 0.4980 
4 80.0 8.0 io 81.5 0.0873 02 5180 
1.14 90.0 9.0 89.7 91.4 0.0873 0.4930 
eso 10.0 1.0 12.9 Pics 0.1745 1.0000 
1.34 20.0 2.0 ZO 7 230 0.1745 0.8140 
1.24 30.0 3.0 29.9 Sor. 2 0.1745 0.6380 
eS 40.0 4.0 39.5 42.9 0.1745 0.5640 
eed 50.0 a. 0 49.3 Be / 0.1745 Os eile 
ee 60.0 6.0 a) el 62.5 0.1745 Oe ort 
eo 70.0 20 69.0 72.4 0.1745 0.5210 
1.14 80.0 6.0 78.9 82.3 0.1745 0.5070 
1.14 90.0 9.0 88.8 92.3 0.1745 0.5000 
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AVG. NO. SHIP SPEED RELATIVE SPEED ALFA PROB. 


LEGS SPEED RATIO UPLEG DOWNLEG DETECTION 
1.40 10.0 Lu Ze PD 39 0.2618 ee OUe0 
3.2 20.0 2 0, ra 24.6. 0.26038 0. SarZ 0 
a6 30.0 3.0 Ze. | 34.0 0.2618 0.6500 
EY 40.0 4.0 3. 0 43.7 0.2618 Oat) 10 
e20 310;3)0 5.0 48.4 oB.5 0.2 Gers 0.5260 
ama 60.0 6.0 Dp ac B3..3 0.26182 0.5270 
iS 70.0 yO 63.1 feo. 2 0.26158 Os ILS 10 
1.14 80.0 8.0 78.0 Be. 2 0.2618 0: 35 0 
LS 90.0 9.0 Gy .9 On3.. 1 0.2628 O58 10 
oo EOE ew Le 16.4 0.3491 ie 000 
en 2 PAB) 10 Zao ho 1 eve 0.3491 0.8600 
a4 30.0 S 0 ZB 34.7 0. 349m 0.6920 
0 40.0 4.0 a7. 8 44,4 0.3491 0.6670 
Lead 5 Ora 9.0 47.5 54.2 0.3491 0.5470 
aS 60.0 6.0 57.4 64.1 0.3491 0.5480 
lee 0 10m 7.0 Gy. 2 74.0 0.3491 O23 290 
1.14 80.0 8.0 wa. 2 ge .9 0.3491 0.5340 
eel 90.0 or0 om. 1 ve. 9 0. 34 0.5070 
1a 10.0 Lee 0 10.7 Fe. 9 0.4363 1.0000 
ees 0 20.0 Zoe) Bre 2 Zr. o 0.4363 0.8570 
iro S10 (0 B70 a. 3 35.4 0.4363 OF 7aa0 
aie 40.0 4.0 oe. 9 45.1 0.4363 0.6430 
eeeles 30.0 eae 46.7 35.0 0.4363 0. SoHo 
Pee 60.0 6.0 56.5 64.9 0.4363 0.5530 
dla 70m O nO 66.4 74.8 0.4363 0.5860 
sO 80.0 oe 0 76.3 84.7 0.4363 0.5360 
ears 2050 oO) 86.3 94.7 0.4363 0.5490 
Lee LOO 10 10.0 i. 3 0.5230 iO C00 
i 0 20.0 2.0 i... 3 2069 05526 0.8890 
R26 30.0 3.0 Zo... 5 S67. 1 0.5226 O. 7200 
1.24 40.0 4.0 36.1 45.8 025256 0.6690 
Wa 1 5070 5) 510, 45.8 oS. / 0.5226 0.6100 
LO G00 6.0 oe. 7 65.6 0:5236 0.5830 
eS 70.0 ee 0 65.6 dae D 0.5236 0. S920 
dive 80.0 3.0 ie. 5 85.4 0.5236 0.5760 
eres 90.0 ese 85.4 95.4 0.5286 0.5540 
ie 7 1Or=0 0 re egies Oo 102 1.0000 
Lge. 20.0 2.0 16.4 Zi 210 065509 0.9010 
1.26 3020 3.0 2.6 S5./ 0. Go 0.8140 
ie2o 40.0 4.0 oe. 2 46.5 0.6109 0.7480 
Leg DOra0 rai0 45.0 36.3 026 Lug 0. 70980 
Lo 60-0 67, 0 54.9 66.2 O56 108 0.6540 
dele? 70.0 ew 64.8 76.2 G6 0s 0.6190 
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MONTE CARLO TECHNIQUE 
PROBABILITY OF RANDOM ENCOUNTER 
NORMALIZED CONVOY ROUTE 


DEFINITE RANGE LAW 

DETECTION RANGE: 30.0 MILES 
LENGTH OF SUB PATROL LINE: 120.0 MILES 
X-DISTANCE TO PATROL LINE: 30.0 MILES 
SUBMARINE SPEED: 10.0 KNOTS 


NO. OF SAMPLES TAKEN FOR EACH SCENARIO: 1000 
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PROBABILITY OF RANDOM ENCOUNTER 
EXTENDED MODEL 
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APPENDIX B 


MAIN COMPUTER PROGRAM FLOW DIAGRAM 


START 


INITIALIZE VARIABLES 
AND COUNTERS 


DETERMINE COORDINATES 
OF LOWER POINT OF SUB 
PATROL LINE 





COMPUTE MAX TIME THAT 
SHIP WILL BE WITHIN 
RADIUS OF DETECTION 

RANGE FROM SUB PATROL 











DETERMINE INITIAL 
SUBSROSITION VECTOR 
FROM SH 





7] 


COMPUTE TIME TO GO FOR 


SUB TO REACH END OF 
PATRCI, FROMALTS 
INITIAL POSITION 


V/ 





COMPUTE RELATIVE DISTANCE 
SUB MOVES ON THIS COURSE 
(LEG) 


COMPUTE FINAL SUB 
POSITION VECTOR FOR 
THIS LEG 














DETERMINE 
CPA FOR THIS 
LEG 


no (detection does 
CPA < R not occur) 


V, V/ 
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yes 


TAKE ANOTHER 





yes (detection occurs) \/ 











Ohi Y ator YOND. R 
FROM PATROL LINE AT 
TIME END OF THIS LEG 
IS REACHED 







no 


10000 legs no 
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yes 


SELECT 
= - NEXT 


no 960 SAMPLES 
TAKEN 
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PRINT RESULTS 
FOR THIS SPEED RATIO 


CROSSING POINT 


AVERAGE NO. LEGS 
SHIP SPEED 
SPEED RATIO 

RELATIVE SPEEDS 
ROUTE ANGLE 
PROB. DETECTION = 


(no. detections) 
(no. samples) 


SHIP SPEED = 90 KNOTS 
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SHIP SPEED < 25 KNOTS 


yes 
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INCREMENT SHIP 
SPEED BY 5 
KNOTS 
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yes 


BASIC MODEL 





INCREMENT a>tan~+(2R/p,) ue 
SHIP SPEED 
BY 10 KTS. 

Vl yes 
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EXTENDED MODEL 


SHIP SPEED > 90 KNOTS 


no 


SHIP SPEED <325 KNOTS 


yes 
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INCREMENT 


SHIP SPEED 
BY 5 KNOTS 
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yes 


EXTENDED MODEL 








CROSSING POINT 
90 MILES FROM 

CENTER OF SUB 

PATROL LINE 






SHIP SPEED> 50 KNOTS 
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> MIEES 


INCREASE CROSS 
POINT DISTANCE 
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20 MILES 
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SUBROUTINE WYEB 
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SUB HEADING "UP" 
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SUBROUTINE CPACHK 


COMPUTE DOT PRODUCT 
(REF VECTOR): REL SBEED VECTOR) 
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CPA LIES BETWEEN 


REFERENCE VECTOR CPA = TERMINAL 
AND TERMINAL VECTOR: VECTOR 
COMPUIES CPA 
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APPENDIX C 


23 IG RIOR IR OR OROR IIE OR a ak ak a ok a 
* NUMERICAL segeaa oe * 


* BA EL * 
SII III OR SiR ROR IR 8 tak 


FOFMAT CONTRCL CCUATER CINITIALIZE) 
SUG SPEED 
V=10. 


SAMPLE SIZE 
JKL=960 


AO. CF ANGLE (ALFA) INCREMENTS 
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GC TC 22 
21 DALF=(1A+3)*5 
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22 U=10. 
CALCULATE POSITION OF BCTTCM CF PATRCL LINE TO 
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YP=XP*TAN(CALF)~(PL/2.}) 
LOOP TO INCREMENT RHC 
DO 801 IB=1,IBI 
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KERR RAK KAKA KK AH RHR HH 
* NUMERICAL INTEGRATION * 
* EXTENDE MOOEL * 
CRRAERAAKKAERK ASKER EREKREEE 
CCMMON/ONE/SGNeYB,ZN 


DIMENSION UV(2) »W1V(2)_,h2V(2) 
IRAVE SO SUWLVEZ) {UW2VED) STEMI ote 


P1=3.1415926536 
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OF PATROL LINE 


TE TRAT THE FIRST (RANDOM) LEG IS TC BE EX= 
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* TO GO FOR SUB TC REACH UPPER LIMIT CF PATROL 


INITIALIZE TOTAL TIME 


MM=] 


OC 


AO 


CGO 


103 


101 


102 


106 


104 


COMPUTE LENGTH OF RANDOM LEG 
GO TO 104 
LOOP TO INCREMENT LEG COUNT 
D 0 IN=1,MM 
CHECK TO DETERMINE INITIAL DIRECTION CF SUB 
IF (SGN) 101,105,102 
INITIAL DIRECTION IS UP$ USE PARAMETERS FOR THE 
UP LEG 
TW=TW2 
CALL COPI (UWV,UR2V) 
CALL COPI(WV,W2V) 
CALL COPI(VV,V2V) 
GO TO 106 
INITIAL DIRECTION IS DOWN: USE PARAMETERS FOR THE 
DOWN LEG 
CALL COPI(UWV,UWLV) 
Tw=Tw1 
W=Wl 
CALL COPI(WV,WLV) 
CALL COPL(VV;V1V) 
REVERSE DIRECTION CF NEXT LEG 
SGN=-SGN 
COMPUTE RELATIVE LENGTH OF THIS LEG 
CALL VMULT(UWV,TEMV, TW) 
DETERMINE VECTOR FOR FINAL POSITION ON THIS LEG 
CALL VADD(R1LV,TEMV,RNV,1) 
CALL CPACHK(RLV, RNV,UWV SRMIN, ING J) 
INCREMENT LEG CCUNT 
SLEG=SLEG+e1. 
DETERMINE IF DETECTION OCCURRED ON THIS LEG 
IF(RMIN-R) 353,11 
NO DETECTION THIS LEG: DETERMINE IF SKIP IS OUT 
OF RANGE BEYOND PATROL LINE 
TN=TN4+T 
IF (TN-TMAX)6 9858 
IF THIS IS NOT THE FIRST (RANDOM) LEG», CONTINUE 
1E(MM—-1)105,9,10 
CONTINUE 
GO TO 8 
LENGTH OF TIME REQUIRED FOR SUB TO TRAVEL ENTIRE 
repiagtAN OF PATROL LINE 


PROVIDE FOR POSSIBILITY OF 10000 LEGS MAXIMUM 
MM=10000 


tial LENGTHS OF SECONO AND SUCCEECING LEGS 


O44 
Orr 
Nm 
it il 
Omt—4 
x 


REMENT DETECTION COUNTER 
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CCMPUTE PRCBABILITY CF DETECTION 
7X JIKL 


P=DET 


CF LEGS FOR THIS SCENARIO 


VERAGE AQ. 


A 
KL 
GENERATE CROSSING CCORDINATE FCR PRINTOUT 


XYP=-YP 


PUTE 
G/XKS 


a UW 


ACH SPEEC RATIO 


RE 
Cowl sW2,XY¥PoP 


PEEC BY 5 KNCTS IF RHC<3, 


SHIF S 
O KACTS 


1 


DE RESULTS FCR TWC 
PAGE 


NCLU 
EACH 


u> | 
< 


IF (I 


E FEEC BETWEEN VALUES OF ALFA 


INSERT COLUMN HEADINGS 


NEW PAGE; 
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=ZTRTWOOWOe WSS Bw 
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LO. STOP 
END 
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* 
WK eK aC a oe ae ae ae ak Oe 


I 
x 


2 IO IOI OO IOI I iO IR OROR i i ag koi ak ak a 
* NUMERICAL INTEGRATICN MOCELS * 
SN. NES 


WK She aie ae oe ae a ae oe a ok ok 


x 


OVOO 


SUBFOLTINE WYEBCXP,YP,PLyRC,BCy,A) 


OOOO 


7 ON PATROL LINE 


Sou 
ECTION OF SUB 


Ow 


OOO 


CCURSE WHEN IT REACHES ENC CF 
)+FL)+PL)9,9,10 


Ee, SUB 
BENE 
Ves ve 


REVERS 
PATROL 
#( 2. *( 


IFCZN 


TC SUB 


LY AT TERMINUS OF PATROL 


TERMINE INITIAL RANGE AND BEARING 
oS 


EVERSE CCURSE CF SUB 


AL(A,RO,B,8C,101) 


Pp 
K 
N 


SUBROUTINE VINS(A,X,Y) 


SERTS THE X ANNO VY VECTCR 
RAY A 


<a 


ION ACL) 


SUBRCUTIANE COPI(A,B) 


NERATES THE VECTOR A SUCH THAT 


N A(1),8(1) 


SUBRCUTINE VAODD(A,BsR,N) 


TINE PERFORMS THE FOLLOWING VECTCR 
I 
A 
A 
R 


OUOOUU 


DIMENSTON A(1),0(1),R(1) 
IFIN) 1,293 


VECTORS 


ADO 
4 T=l,2 
D=A(TPOB(I) 
URN 


I 
t 


Or" is 


a) 
4&4 R 
R 


TRACT VECTCRS 


2 
S=B(r) 


CH «tone 


SU 
j= 
A( 
N 


S 
I)= 
TUR 


DC 
R ( 
RE 


1 
5 


TE DCT PRCOUCT CF VECTORS 
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SUBROUTINE VMULT(A,B,X) 


THE 


BY 


THE VECTOR A 


A(1),8(01) 


SUBROUTINE VVAL(A,VL»,UV,ANGLE »M) 


CCW FROM X-AXIS 


L TO A 
RED 


UPPLIED BY THE CALLING PROGRAM) 
fete 
ASU 


INE DETERMINES THE FOLLOWING 


OVOVVNIO 


A( lL) ,UV(1) 
2e3° 
5 ,6,7 


F VECTOR ONLY 


C 
1 


AND ANGLE OF VECTOR 


ce GENERATE UNIT VECTOR 


oF 


CPACKK(R1,RNZUWV,RCPA, ING J) 


SUBRCUTINE 


OOO 


Owe O 
aio ez 
Ui bjO 

Qa Vi UL 
Lick oes 
Lo OU Q 
lirik We 
re i qd ct 
Wthraiwn 
ate SL Il 
we Oo es OC 


OUOOOO 


eee FIC LI »RN(1) pUWV(1),X02),Y(2),CPA(2),2(2), 


CA/TWO/Y,N 


DIME 
1FT(2 
CCMM 


R1*UWV=X 


UCT: 


sel CN THIS LEG; CCMFUTE TRE COT FROD- 


CC 
*UWV 
UbV 


ot 


“ 
( 
’ 


PA WIL 
CT: RN 
ADD(KN 
Pele sveu tes 


T CCCUR CN THIS LEG: CPA=R1 
A,yX»C,100) 


=RN 


CPA 


N THIS LEG; 


Oe 


X*¥UWV=Z 


ASSUME LEG HAS LENGTH < 250 MILES 
98 


BY Ol MILE AND FIND CPA 


RECECIANG VECTOR 


Y »P-2.) 


X,1) 
CHECK LAST 2 MILES CF ORM IN O/f1-MILE INCREMENTS 


T=l,21 


P=0.1%*] 


CT PRCCUCT 


p 
9 
9 


ENT 
UWV»Y,?P) 
C REFERENCE VECTOR 


D0 T 
DDO(R1»Y¥,X,1) 


MENT CRM BY 2 MILES 
(UW 

COMFUT 

VACO(X 


MPUTE D 
OD(X,UWV,Z 40) 
INCREMENT CRM BY 0.1 MILE 


USE 2-MILE EANCREMENTS 


Zot 
AD 
VA 


P= 
CALL 
CALL 
DC 8 


RROR_ EXISTS IN PROGRAM IF NEXT INSTRUCTION IS 
XECUTED 


E 
E 


OO 


AMPLE*,1551X_,*GENER', 


IF NEXT INSTRUCTION IS 
Xe°'CF SAMPLE*,15,1X,°IS 
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250 MILES 


“y 


oO 


TIE 


APPENDIX D 


He WK ee fe ok aie ak a ake oe ic ok ak oe oe ooo ok 
* GANDOM SAMPLING * 
* MIDPOINT CRISSING * 
We nde fe ie a he ok ae ee ak ak ae te aioe ie ae a 


COMMON/OBR F/SSGN,YB 


NIMENSTCE! UIVO2) pWLVC2)gW2VO2) pVIVEI2Z) ep V2VE2)9R1VI2)9 
LEMNVO2 Dg UNTV( 2), UW2V(2) 6 TEMVI2 Sg WV C2), VV02),UWNY(2) 
IMITICLET7E URN GENERATOR 


POST THUR (C) 
mae. 1 TSS ss 
SUR SPEEN 
veut se 

SAMPLE SI7E 
JKt=1ecP 

Mt. GF FNGLE WALFEA) INCREMENTS 
Ret =12 

NO, OF FHO INCREMENTS 
TRI =O 


CERNE Oath SURBVAZCTOCRS 


CF WAIN IGN 5 ice V.) 
CALL VINS(V2V,0.9-V) 
SIZ] TF DETECT IGNAPADTIUS (DEFINITE RANGE LAW} 
F=3C. 
LENGTH TF SUR PATROL 
PL=120. 
TNITEAL X-PNSTTION NF PATROL LINE 
XE=R 


WEI TE CSE COIR, OL XP a Va JKL 
INITIAL ANGLE 
AL E=C,A 
LOCP TO INCREMENT ALFA 
[rr PNM [fL=l,ITAl 
INITI“ZL SHIP SPEED 
1l=10, 


CebGuibnie POSITION Ch AROT TOMS OF DATRE LINE 70) 
INSURE CROSSING AT CENTER OF PATRNL 


VF=HXPRTENCALFI-CPL/2. ) 


100 


CAO 


oy 


ANAOHAAA AQAOQAAAGY 


OOCE 


LOOP TO INCREMENT RHO 
DO 801 IB=1,IBI 

INITIALIZE LEG COUNT 
SLEG=0. 


COSA=COStUALF) 
SINA=SINCALF) 


DEFINE VECTOR OF SHIP VELOCITY 
CALL VINS(UV,U*COSA,U*SINA) 

COMPUTE SHIP/SUB SPEED RATIO 
RHO=U/V : 


INITIALIZE SUB DIRECTION FOR DETERMINISTIC 
CALCULATI : 


SGN=1. 

INITIALIZE DETECTION COUNT 
DET=0. 

COMPUTE MAXIMUM TIME FOR CALCULATIONS 
TMAX=(XP4R)/(U*COSA) 41. 

DETERMINE RELATIVE SPEED VECTORS AND VALUES OF 


LIZE 
ATION 


MAGNITUDES AND ANGLES 
CALL VADD(V1V,UVeW1Ve~1) 
CALL VADD(V2V,UV,W2V,-1) 
CALL VVAL(CWLV,W1,UW1V,BW1,111) 
CALL VVAL(W2V,W2,UW2V,BW2 9111) 
DETERMINE IF DETECTION OCCURS 
EACH COUNT OF THE LOOP INDEX REPRESENTS ONE SAMPLE 
DO 8 J=lyJKL 
DETERMINE SUB POSIT (UNIFORMLY DISTRIBUTED) 
XP: X=-DISTANCE TO PATROL LINE 
YP: Y-COORDINATE OF LOWER LIMIT OF PATROL 
PL: LENGTH OF PATROL LINE 
OUTPUTS 
RO: RANGE TO SUB*S INITIAL POSIT 
BO: BEARING OF SUB*S INITIAL POSIT 
YB (VIA COMMON/ONE/): INITIAL Y-COORDINATE 
Riv: INITIAL SUB POSITION VECTOR 


CALL WYEB(XP,YP,PL,RO,B80,R1V) 
TIME TO END OF FIRST (RANDOM) LEG 
IF (SGN) 700,105,702 


TIME TO GO FOR 
IF INITIAL SUB 


T1=(YB-YP)/V 
GO TO 703 


UB TO T QF PATROL LINE 
ELOCI = 


S 
VELO 


TIME TO GO FOR SUB TO REACH UPPER LIMIT OF PATROL 


101 


Kee 


Lea 


LIF 


TIS(YP+EL “VAI /Y 
INITIALIZE TOTAL TYME 
Th=C, 
INDIC¢TE THAT THE FIRST (RANDOM) LEG IS TO RF 
“XAMENES 
MMS] 
COMPUTE LENGTH OF RANDOM LEG 
eC Tr yA4 


LCEP TE SINCREMENT LEG CGBNT 
Oe TONS 

CRECK Tl DETERMINE INITIAL DIRECTION OF SUS 
TEC SGN TCI, 195.102 

[NT DEAL PIRECTIONSITS VPs USE PARAMETER Sage rai 

UF LEG 
W2 
COFT CUWV,UnW2V) 
COPTEUWV,W2V) 
CONFT(VV,V2V) 
196 


po a TS NOWNs USE PARAMETERS FOR THE 


DONE D4 
JPY 1 


n 
n 
1 


eW1V) 
»V1V) 
SeGrran Cee NEXT Les 

COMPUTE RELATIVE LENGTH OF THIS LEG 
CALL VMULT(UWV,TEMV,TW) 


CE TERMPNe VECTOR PER FINAL POSITION ON THIS Leo 


CALL VADDUFLV,TEMV,RNV, 1) 
CELL CPACKK(RIV,&NV,UWV sRMIN, ING J) 


INCE SMERT LEG oun 
SLEG=SLEG+4+1. 
DETTE MINE Tr VemecCt ren UCCURS EG ONG Ss LEG 


SS het ae s3ell 
TheTh 

TEL TNOTMAX) 4 848 
TF(MM=1910&,9,10 
CONTIN 

Cot B 

T=OL/V 

MM=}TOOCO 


ik0!2 


0 

2 

? 

9 
DETECTION! /T5, *LEGS 
LLY LIKELY) OF SUB 


ROL LINE 


EG,U,yRHCyW1] pW2,ALFyP 
B 
W 
0 
5 
N 
Q 
F 


=a = Mme OMRZU Win 
Kee MN BRePMree es © cai. 
e aI™G OC mete ™ eg ZWWY e 
Om ¥OO OR GATusOelwwWwo> + 
FANO¢WwwI7W @ en ¢ Se SEO Ol 
SSK POM MYO eOwUlw OWA WW 
He WZM™”N~Cet wk WNDFeE afb 
me OO Ih HW def ele | AQ. < 
Hh be HO +e =H ROMNMw MOAWWNS AO. 
FN &erSoMOW HS DHOULAE EW Wea OO 
SBOWO NINH WAOdON | OQewrvnen OFZ 
me OOUKaAxKEDPDBUCUN]| of eee LWW 
MNO Siw Of © 
+ (mc st NOOO m= 
© Cc Ooo oOo 
oo wm ODO 0 = 


Y 
R 
E 
S 
p 
E 
D 
U 
2 
SUBROUTINE WYEB(XP,YP,PLyRO, BO, A) 
I 
0 
L 


WOOOUOO 


DETERMINE INITIAL RANGE AND BEARING TO SUB 
103 


DETERMINE SUB POSIT RANDOMLY 


YB=Y¥P+URN(1) *PL 
DETERMINE INITIAL DIRECTION 
DOWN 


UP 
7 SGN=-l. 


cab ty aay 


yay 


— 
2 


x» 


— 


SURFCUTINE COACKHK(RLILZRN,UWV,RCPA,IN, J) 


FIMERSTAM RLCLIVARNOD DP ,UWV(L)SX02),¥(2),CPA(2),72(2) 
COMMCR/TWO/Y AN 


ky=REFERENCE VECTOP 
RM=TCEMYINAL VECTOR 
UWV=UNYIT VECTOR CF DRM 
COMFUTE NOT PRODUCT: RL *UWV=x 
CLLL VLAN( RL eUWV,X-D) 
PEPCK O1L) Tia 2a 
CPA WILt OCCUR ON THIS LEG: COMPI'ITE THE ADT 
PFOOUCT: RN*UWV=X 
CfLt VADDIRN,UWV,X,9) 
TFCXC LIV 3» 304 
hate WILt NOT OCCUR ON THIS LEG: CPA=R] CAMOISTE 
Fs | 
CELL VVSLIF1L,RCPAX,0,190) 
RE TUF A 
aes WAS NOT REACHEG ON THIS LEG: CPA=RN: COMPLITE 
CALL VVALCFEN,ZROCPA,X,0,109) 
QFTURN 


NET PM) Nc CR eeASSUMETLIEEGSHAS LENGTH< 2504ILES 
DSS em lege) 
USFeHS Ce S4rCe INCREMENTS 
P=1/2. 
IMCPEMENT DPM OY 1/2 MILE 
CAttL VMULT(UWV,Y,P) 
AC'D TO BERERENCE VFEROR 
CELL VADO( EL» Y¥,X%_1) 
COMPUTE DOT PRNDUCT xX NOT UWV=Z 
CALL VADER (X,UWV,7,0) 
1E @ CHPNGES rN BECOMES ZERO VECTOR X=CPA: 


OQTHFFWIS7 TRY 
LE OZ” .BgSegh 
CONTINUE 

DPM LFG>259 MILES 
WPITE(6,6C) IN, J 
FOPMAT(S/ICX,*LEG* » 1X pl&e1Xe'OF SAMPLE*,15,1X,'°IS 
GREATER THEN 259 MELEStS7/) 
RE TURN 
CALL VVAL(X,RCPA,Y,0,100) 
PF TURAN 
ENN 


104 


AOIOOOASY 


uw 


ON 


SUBROUTINE VINS(A,X,Y) 
THIS SUBROUTINE Se THE X AND Y VECTOR 


COMPONENTS INTO ARRAY 
DIMENSION A(1) 
A(1)=X 
A(2)=Y 
RETURN 
END 


SUBROUTINE COPI(A,B) 


THIS SUBROUTINE GENERATES THE VECTOR A SUCH THAT 
FOR A GIVEN VECTOR By, A=B 


DIMENSION A(1),B(L) 
D0 1 1=1,2 
A(1)=B(15 

RETURN 

END 


SUBROUTINE VADD(A,B,R,N) 
aS weno DOES THE FOLLOWING VECTOR 


am = (flee 00 00 


ADD VECTORS 
4 1=1,2 
=A(T$4+B(1) 

RN 

SUBTRACT VECTORS 


1=1,2 
=A(1}-B(1) 
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Cet) ah 0d 


G 


(N 


SUPPCUTINEG VMULT(C4,B,X) 


THIS SUPR OUT INET IP 
SC4LELR XY 7) FORME THE ¥ 


DIMENCTON ACL), BOL) 
ae yO 

ac py=xeelT) 

Cr TUPN 

CAD 


TES THE VSCT OP eee AE 
CTAR RB 


amma 


SUPFECUTINEG VVALCA,VL»UV,ANGL= 9M) 


THIS SURROUTINE DETERMINES THE FOLLOWING 
A=TNPUT Vearer (SUPPLIED RY THE CALLING PROGRAM) 
VLEMAGNITUDE OF A 
UV=URTT VECTOR PARALLEL TN A 
ENGLO=ANGLE OF A MEASURED CCW FROM X-AXTS 


CAME MSs me ay) 
Let Mara, 
TEUMoLIC Ee tAe eT 


COMPUTE VECTOR LENGTH ONLY 
=$ 


VLHSORTC AC) )**A4A(2)*H2) 
PF TURN 


COMPUTE SNGIE OCR YECT@R ONLY 


ANGLE=ATAN2Z(A(2),A(1)) 
RETUPN 


COMPUTE BOTH LENGTH AND ANGLE QF VECTAR 


VE=SORTCACLI##24¢A(2 142) 
Gr TN 


GENERATE UNIT VECTFER ONLY 


COMOUTF ANGLE AND GENERATE UNIT VECTOR 
F= mae Vaeeceniae? 
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AOOQO 


KKAKKRHRRARRRAARARRAKAAKRR RK ARE 
* RANDOM SAMPLING x 


* NGRMALIZED CONVOY ROUTE * 
WR RR RR RR RRR RK OK BK EK HER RK 


COMMON/ONE/SGN,YB 
DIMENSION UV( 2) 4WLV(2) oW2V(2) pVIV(2)9V2V(2).R1VI2), 
(2), TEMV(2$,WV(2)5VV(2),UNV(2) 


LIRNV(2),UW1V(2),UW2V 


REAL NORM 

INITIALIZE URN GENERATOR 
POSIT T=URN(0) 
PT=3.1415926536 

SUB SPEED 
V=10. 

SAMPLE SIZE 
JKL=1000 

NO. OF ANGLE (ALFA) INCREMENTS 
TAI=17 

NO. OF RHO INCREMENTS 
IBI=9 

DEFINE BOTH SUB VECTORS 


CALL VINS(V1V,0.,V) 
CALL VINS(V2V,0. ,-V) 


SIZE OF DETECTION RADIUS (DEFINITE RANGE LAW) 
R=30. 
LENGTH CF SUB PATROL 
PL=120. 
INITIAL X-POSITION OF PATROL LINE 
XP=R 
WRITE (6,600)R_ PL XP eVeJKL 
INITIAL ANGLE 
ALF=0.0 
LOOP TO INCREMENT ALFA 
DO 800 IA=1,IAI 
INITIAL SHIP SPEEDO 
U=10. 
LOOP TO INCREMENT RHO 
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OID 


ee) 


OY 


cy 


AOAC Velo 


YI a 


Cr 90] gipeiae ou 
INT TEAC Teer GUN T 
SoG —or 


(@<pA=COS( se) 
STNA=SINCALFEY 


PEFINE VECTOR OF SHIP VELOCITY 
CALl VINSCUV,Z,UeCOSA,U*SINAD 

COMBINE SHI P/SURSS Eee Amr) 
REO=HUSV 


INTITIALIT2 
CAC Ue ar at 


SGN=1. 

INITIALIZE CETECTICN CCUNT 
Mt T=C, 

COMPUTE MAXIMUM TIME FOR CALCULATIONS 
TMAK=(XP#PI/S(UECNSAD 41. 


Ce TSR ete EN SPEED VECTORS AND VALUES NF 
MAGNITUPES AND ANGLES 


canes NIRECTION FOR DETERMINISTIC 


CALL VADN(VLV,UV»WLV el) 
CALL VADDCV2V,UV»W2V,-1 ) 
CALL VVALOWLV, WL UWLV,8W1 lL) 
CALL VVALCW2VyW2,UW2V,BW2,111) 
NOPMALT7& CONVOY TRACK FOR EACH UNIFORMLY 
DISTRIBUTED SUBMARINE STARTING POSITION 
SAMPLE SI7E FOR NORMAL IS CETERMINEN BY THIS Lane 
np? Met, 1005 
YE=NORM(26, )—(PL/2.) 


( 
DETERMINE TF DETECTION OCCURS 
EACH CCUNT CF THE LOOP INDEX REPRESENTS INE SAMPLE 


Cae J=1,J 
DETERMINE SUR POSIT (UNIFORMLY DISTRIBUTED) 
Xe: ~DISTANCE TO PATROL LINE 
YE: v=COORDINATE OF LOWER LIMIT CF PATROL 
PL: LENGTH OF PATROL LINE 
OUTPUTS 
FOS RANGE TO SUB'S INITIAL POSIT 
BC: BEARING OF SUR'S INITIAL POSIT 
VR (vrs COMMON/ONFE/): INITIAL Y-CIORNINATE 
PIV: INITIAL SUB POSITION VECTOR 


CALL WYERCXP,YP,PL,RO,BN,RIV) 
TIMF TO ENA SEF FIRST (RANDOM) LEG 
TF (SGN) 7€C »105,702 


TEM VT GsG oO eGR 


SUB _ TO RF QF PATROL LINE 
IF INITYAL SUB VELOCITY 
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CIc7 


Cy 


a 


i) 


oy 


ca 


2 


3 


&¢ 


Y 

i. 
TIME TC GO FER SUR 1G PEACH UPPER ETM ieee snr Ot 
ENS 

Y 


TL=CYP+PL—-Y8)/V 

T=7T) 
IPYITIALITZE TCTAL TIME 

= 
INMICLTS THAT THE FIRST (RANDOM) LFG ITS TO RE 
FXAMINER 

MM=]} 


COMPUTE LENGTH OF RANDOM LEG 
GE TO 1¢4 

tCGP TP ENCREMENT LEG COUNT 
nr 3c TMN=],MM 

CHECK TO NETERMINE INITIAL OIRECTION OF SUB 
VE USGA tl oa oe 


eres PLOT CT HGN 1S UP. WS] PARAMETERS Sale te 
ae ce 


DAVE D 4 
Atl if» 


TTON RE NE HG 
SGN=-S6N 
COMSUTE Sst Vee eon a eet Ss LEG 
CALL VMUE TOUIWV, TEMV, Th) 
Pete reine VEC R ECR sein eee esl 1 1TON GON Silt Sere G 


FALL VADP (FLV, TEMV,RNV,1) 
CALL CPECHK(RIVZRNV,UWV,RMIN, IN, J) 


INGE EMERY EEC er Guns 
SUEG=HSLFEG+4+1. 

PUT PRIN T= DETECTION eV ERED INN TH LS she G 
T FOES ate e931! 
ek 
IF(TN MAX 1698 9A 


02 


EGsUs,RHC Wl »W2,ALF,P 


8 


MM=TCOCC 
162 
cT+l. 


Tew] 
29 CONTINUF 


Tew? 


= 
=> 
— 
=> 
= 
‘ 


FUN) TCS ean 


] 
CONTINUE 
CC Ge 
T=PL/V 
TW2 


> 
9 
194 Th 


l 
oS 


a mle MG 
* Us TC) oe 

ee ZOOL 
OS OT Ulic 
eit oe Wed 
eR o. WO. ua 
~~ ad ey! 

entice fC a 


CON & | lCe 
CAN 0 OOK TU 


DemeECTION'/T5,°LEGS 


DCWNLEG*//) 
o2 e2XK oe BCL K FS Dg 2K FS ely 3X ap F be 4s Xy F6L4) 


PRAB,. 


ALFA 


IN UPLEG 


oN we Le CU LN CG. Tv,keu — 
ie eee 
Wiest atcl ele yu, 
Cre e+e FNS Ce]! MOwWWMZu 
Tes gee 

aC wo S 
1 Se 


dwt 
aA < 


eLumeao ACC 
Fe yr eee CFs 
wh ee OU 


amNY & 


Qi: 
CC) Cae 
Shien Gra! hae 


ir 
CG: 


SC ane 


SURROUTINE WYEB(XP,YP,PL,RO,B0,A) 


PCY LIKECNI ae Sus 


ROL LINE 


es 


PA 
EB 
Oy) 


J 
: 


N 
8 
( 


UOUGO 


DETERMINE SUR POSIT RANDOMLY 


YP+UFN( 1) *PL 


{2y 


INT TIAL RANGF AND BEARING TN SUB 
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—™ 
- 


AL(A,RO,R,BO,191) 


DETFRMIRS INITIAL DIRECTION 
XF 
YR 
VV 


4 TFCURN(1)-€.5)5596,7 
é G6 TO 4 

DOWN 

DE TFRMIR 


UP 
7 SGN=-)]. 


5 SCMN=1. 


GO Te 8 


YR 


OAANY 
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SUBROUTINE CPACHK(R1,RN,»UWV,RCPA,IN,J) 


DIMENSION R1(1),RNO1) SUWVO1L) XC 20,V02),CPA(2),2(2) 
COMMON/TWO/SY,N 


RI=REFERENCE VECTOR 
N=TERMINAL VECTOR 
UWV=UNIT VECTOR CF DRM 


COMPUTE DOT PRODUCT: R1*UWV=X 
Tio See 


9 


L OCCUR ON THIS LEG: COMPUTE THE DOT 
> RN*UWV=X 


O 
~F Cx = 
— WS Om re 


9 
L 
ir 
(RN»UWV,X,0) 

9 4 

LL NOT OCCUR ON THIS LEG: CPA=R1 COMPUTE 
VAL(F1,RCPA,X,D,1090) 

WAS NOT REACHED ON THIS LEG: CPA=RN: COMPUTE 


VAL(URN,RCPA,X,D,100) 


Zc pp Fe pp ewes OP we 


OETERMINE CPA: ASSUME LEG HAS LENGTH< 250 MILES 
DO 5 1=1,500 
USE HALF-MILE INCREMENTS 
P=I1/2-. 
INCREMENT DRM BY 1/2 MILE 
CALL VMULT(UWV,Y,P) 
ADD TO REFERENCE VECTOR 
CALL VADD(R1,Y,X,1) 
COMPUTE DOT PRODUCT X DOT UWV=2 
CALL VADD(X ,UWV,Z,0) 


IF Z CHANGES SIGN, OR BECOMES ZERO VECTOR X=CPA: 
OTHERWISE TRY AGAIN 


IF(Z(1)) 59656 
CONTINU UE 


DRM LEG>250 MILES 


RITE(6,60)IN,J 
ORMAT(//10X, TLEG' 4 1X914,1X, "OF SAMPLE, I15,1X,°IS 


1a bis 


AOOODAAY 


P THEN Zo FILES * 77) 
VAL (X,RCPA,V,D,109) 


oLaG 
URN 
L V 
UR 


N 


MIADD 
2nmPND 
2a 


SUBFRCUTINE VMULT(4A,B,X) 


THIS SURROUTINE MULTIPLIES 
SCALAR X TO FORM THE VECTOR 


N A(1),801) 


pate VECTGR A ®Y THE 


mt NN) 


SURBPCUTINE VVAL(A,VL,tIV,ANGLE »M) 


VALUES ne AE DETERMINES THE FOLLOWING 

=INPUT VECTOR 6 celal BY THE CALLING PROGRAM) 
VL=MAGNI TU DE OF 

UV2UNIT VECTOR PARALLEL 0 


TO A 
ANGLE=ANGLE OF A MEASURED CCW FROM X=-AXIS 


DIMENSTON aren. ee 
IF (M-100)2, 
TE(MoL1O}8 (607 

COMPUTE VECTOR LENGTH ONLY 
Geet ont a Pn CA Fikes 
RETURN 

COMPUTE ANGLE OF VECTOR ONLY 
aE cae LE= yee ene eI 


COMPUTE BOTH LENGTH AND ANGLE OF VECTOR 


VL=SORTCA( 1) **®24A(2)**2) 
GO TO 2 


GENERATE UNIT VECTOR ONLY 
Vi=SORT(A(1)##2¢4( 21882) 

UVC CT yeati/ve 
RET 


CCMPUTE ANGLE AND GENERATE UNIT VECTOR 
ANGLE=ATAN2(A(2),A(1)) 
GG TC 6 
END 
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SURROUTINE VINS(AsXsY) 


THIS SUBROUTINE eats THE X AND Y¥ VECTOR 
COMPONENTS INTO ARRAY 


nee SION A1) 
) 


N 
UR 


I 
( 
( 
E 
N 


DIM 
Atl 
A(2 . 
RETURN 
END 


SUBROUTINE COPI(A,B) 


THIS SUBROUTINE oe peace THE VECTOR A SUCH THAT 
FOR A GIVEN VECTOR By A= 


MENSION A(1),B8(C1) 
, =1,2 


I 
C = 
(1y=BCTS 
E TURN 

ND 


MVYPagod 


SUBROUTINE VADD(A,BeR,N) 
ape SUBROUTINE DOES THE FOLLOWING VECTOR 


PERATIONS: 
TF NO: R=A+B 
IF N<O: R=A—B 
IF N=0: THE DOT PRODUCT OF A AND B IS FORMED AND 
RESULT IS PLACED IN THE FIRST MEMBER OF ARRAY R 
DIMENSION A(1),B(1),R¢1) 
IF (NJ 12293 
ADD VECTORS 
DO 4 I=1l,2 
R(YP=A(T$4B(1) 
RETURN 


SUBTRACT VECTORS 


5 [=l,2 
=A(T$-B(1) 
URN 
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APPENDIX E 

A. GEOMETRIC INTEGRATION 

Figure 7 is a graph of a family of curves plotted from 
Table E-la which represents the conditional probability 
that detection occurs as a function of the crossing point 
displacement from the patrol line center. The normal 
density function, with mean zero and standard deviation 20 
miles, is also plotted in the figure. The upper set of 
numbers along the abscissa represents the normalized values, 


ee: | 
x' = x/20 


where x' is the normalized value of x. The product of the 
normal density curve with each of the conditional proba- 
bility curves is tabulated in Table E-lb, and plotted as 

a family of curves in Figure 9. If a grid of 0.1 milli- 
meter squares is superimposed on this figure, the area 
represents one-half of the integral of the function. The 
area under each curve was determined by the grid-counting 
method, and was also determined by use of the Amsler Polar 
Planimeter. Each curve was traced five times with the 
planimeter stylus, and an average planimeter dial reading 
taken. The results of these two methods of geometric 
integration are tabulated in Table 9. The average plani- 
meter dial readings were normalized to 0.998 (the value of 


P(D) for p = 0) using the following equation: 


LT 


TABER E—. 
CONDITIONAL PROBABILITY OF DETECTION 


GIVEN CROSSING POINT IS DISPLACED 
Y MILES FROM PATROL CENTER 


a. NO ASSUMBTION ON DISTRIBUTION OF CROSSING POINTS 


oD 79 yp S4 oe S635, FOG 869885 SSARP OJ 3 38E 
omen cli TL I6O™ sO 2 slays ro) eeincion 


ED Yp = y) 

4 SHIP-SUBMARINE SPEED RATIO (0p) 

0=0 o=1 o=1.5 O=2 o=3 9=5 o=9 0 =00 
0 1 i 9354 .8000 .6250 .5104 .5021 .S@bE 
5 1 1 -9200 ./7844  .62590 .35260 “502. Peace 
10 1 1 -8687 .7635 .6177 .5469 .502]7 25008 
nS a -9781 .8062 .7125 .6104 .5604 .5292 s5000 
20 1 9365 .7437 .6500 .5969 .5604 .5292 .,5006 
ZO 1 -6844 .6812 245675 §<5552 «5469 .5292 =s0c3 
30 1 so2lo 4.64187 25302 .59985. 25052... .5010 Sno 

1 

ik 


b. NORMAL DISTRIBUTION (STANDARD DEVIATION 20 MILES) 
OF CROSSING POINTS ASSUMED 


P(D|¥,, = y) 
* 
= SHIP SUBMARINE SPEED RATIO (p) 
eo o=1 o=1.5 O=2 o=3 o9=5 o=9 Q=00 


OD eogeo. ~seieo..5/31 3891 2.24092. hegeec. 2003 o.oo 
> -«386/ .386/7 .3681 .3033 «<7 op203@— .1942—— wioae 
JOP ooc «=. ool SSO o9 gees 6S) CCG. ey oo lei 
15 .301Lls 2.294% .242:85 me2l4oe wAlS38 2diG87ee01574 .1506 
20 .2420" Va2zzeo. FESO = 20a73. wiead. ..Jj500 «1281. ize 
25 «1826 6095 oolet44 = 1.1073 «8 el Ode. 09909 0966 ae OOS 
30 BIi2y> .-L064 “0812 0687 20665 “Sess 20649  20Gae 
50g) .Gi:75- 020 Ome0067 20064 S006 lie. 00S 9005 Sem. 0058 
70 #0009 .000G 0002 .000Z2 S004 Sate? .0002 sOuUs 


*Extracted from CRC Standard Mathematical Tables, 16th ed., 
19637, pp. 549-555. 
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Area under Curve = (0.9980 x average reading)/762 (2) 


where the constant 762 is the average planimeter dial 
reading obtained by tracing the curve for 0 = Q. 
The values contained in Table E-2 for the column under the 
heading "Planimeter P(D)" were computed by equation (2) 


using a Friden electronic calculator. 


TABLE E=2 


RESULTS OF TWO GEOMETRICAL 
INTEGRATION METHODS 


ae Lae P (DETECTION) _ 
| GRID COUNTING AMSLER POLAR 


PLANIMETER * 





0 0.9980 0.9980 
1 0.9415 0.9319 
1.5 0.7895 0.7922 
0.6816 0.6828 
3 0.5781 0.5797 
5 0.5233 0.5234 
9 0.5008 0.5011 
© 0.4885 0.4883 


*The planimeter dial readings were normalized so that 
762 increments on the dial represented the known area 0.9980. 


1 
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APPENDIX F 


LIMITING CASE FOR P (DETECTION) 


eves 


uo 
R= 30. 
Eae— ae ie 
YP=40. 
WRITE (6 ,600)R,PL 
Z=YP+PL 
DO eas, 41 
IF (YP-R)3,2,2 
2 P=0 
GO Te) 4 
3 IF(YP+R)5,5,6 
5 P=(2.*R) /PL 
GO" 10 4 
6 Y=AMIN1 (ABS (YP) ,ABS(Z) ) 
S=(YP*zZ) 
W=S/ABS (S) 
P= (R-W*Y) /PL 
4 J=J+1 
IF (J-6)7,8,9 
0 
WRITE (6,602) 
WRITE (6,601) YP,P 
YP=YP-2.5 
STOP 
FORMAT (1H1,15(/) ,20X,' LIMITING CASE FOR P (DETECTION) '/ 
1//5x%, "DETECTION RADIUS (DEREINITE RANGE LAW): ',F430775 
2X, ~' PAPROLeY LENGTH go” ,F4..07 /77X,'YP' ,8X,'P(D)'//) 
601 FORMAT(5X,F5.12,/5X,F6..4) 
602 FORMAT (/) 
END 


Oo 


OO tf ~J 
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Ae 1 


LIMITING CASE FOR P (DETECTION) 


DETECTION RADIUS (DEFINITE RANGE LAW): 30. 
PATROL LENGTH: 120. 


YE P (D) 
40.0 Oro 
Seo O70 
3570 0.0 
S25 O70 
3Onma0 O20 
Leet 0.0208 
Paap, 0.0417 
oc 0.0625 
20.0 OZ 0625 
eS 0.1042 
oO O2i2Z50 
i255 0.1458 
10:0 O2N6GH 

7.5 0.1875 

oa O20 ee 

Zao 0.2292 

0.0 Oza 00 
ei S 0.27038 
=506 0 Onc 
hie SD Oe ks 

—1O.30 022555 
= Zero 0.3542 
ane O27 50 
Sli 023956 
S200 0.4167 
=e 0.4375 
2. 0 0.4583 
Seis D 0.4792 
SoUR0 Oe 000 
eo uae) 0.5000 
OOO 0.5000 
=O eS 0.5000 
-40.0 0.5000 
-42.5 0.5000 
=45..0 0.5000 
—47.5 O3000 
= 30150 0.5000 
=O 0.5000 
=O 0.5000 
=e 0.5000 
=6.0'.0 0.5000 


iro 


Uo Ship-Submarine Speed Ratio: © 
Submarine Patrol Line Length:120 

0.8 miles 
Detection Range (Definite Range 


Law): 30 miles 





0 NEO 20 30 40 50 60 70 80 90 
Crossing Position in Miles (from Patrol Line Center) 


Figure F-l. Probability of Detection vs. Crossing 
Position (Perpendicular Tracks) for the 
Limiting Case 
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